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ABSTRACT 

We present an analysis of UBVRI data from the Selected Area SA 141. By applying 
recalibrated methods of measuring ultraviolet excess (UVX), we approximate abun- 
dances and absolute magnitudes for 368 stars over 1.3 square degrees out to distances 
over 10 kpc. With the density distribution constrained from our previous photometric 
parallax investigations and with sufhcient accounting for the metallicity bias in the 
UVX method, we are able to compare the vertical abundance distribution to those 
measured in previous studies. We find that the abundance distribution has an un- 
derlying uniform component consistent with previous spectroscopic results that posit 
a monometallic thick disk and halo with abundances of [Fe/H] = —0.8 and —1.4, 
respectively. However, there are a number of outlying data points that may indicate 
contamination by more metal-rich halo streams. The absence of vertical abundance 
gradients in the Galactic stellar populations and the possible presence of interloping 
halo streams would be consistent with expectations from merger models of Galaxy 
formation. We find that our UVX method has limited sensitivity in exploring the 
metallicity distribution of the distant Galactic halo, owing to the poor constraint on 
the UBV properties of very metal-poor stars. The derivation of metallicities from 
broadband UBV photometry remains fundamentally sound for the exploration of the 
halo but is in need of both improved calibration and superior data. 
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1 INTRODUCTION 

The three-dimensional abundance distribution of the Milky 
Way field stars is a key constraint on models of Galactic for- 
mation (see, e.g., the reviews in Majewski et al. 1993; Free- 
man & Bland-Hawthorn 2002). The observed abundance dis- 
tribution is the convolution of the enrichment history of the 
Milky Way and the dynamical evolution of its stellar popu- 
lations (as reflected in their spatial and kinematical distri- 
butions). Specific models of galaxy formation make speciflc 
predictions about the abundance distribution measured in 
situ from the fleld stars that can be tested against large 
catalogues of abundance measures. For example, the hy- 
pothesis that the Galactic halo was formed in a global col- 
lapse was primarily supported by stars on more radial orbits 
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(and therefore reaching higher z heights above the Galactic 
midplane) being more metal-poor than stars on planar or- 
bits (Eggen, Lynden-Bell & Sandage 1962). However, the 
global collapse theory was unable to account for the lack 
of an abundance gradient in the Milky Way globular clus- 
ters (Searle & Zinn 1978), leading to the hypothesis that 
the Galaxy formed at least partly through the merging of 
larger systems - a hypothesis now well-supported by stud- 
ies of stellar kinematics and spatial distributions (Majewski, 
Munn & Hawley 1996; Ivezic et al. 2000; Yanny et al. 2000; 
Vivas et al. 2001; Siegel et al. 2002, hereafter S02; New- 
berg et al. 2002; Yanny et al. 2003; Majewski et al. 2003; 
Rocha-Pinto et al. 2004; Belokurov et al. 2006; Grillmair et 
al. 2006; Vivas & Zinn 2006; Juric et al. 2008) as well as the 
properties of globular clusters (see, e.g., Marin- Franch et al. 
2008). 

The abundance distribution is of particular importance 
to understanding the nature of the Galactic thick disk. An 
abundance gradient in the thick disk would favor a scenario 
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in which the thick disk was formed either in the slow late 
stages of the early Galactic collapse or the gradual kinemati- 
cal diffusion of disk stars. On the other hand, a monometallic 
or irregular metallicity distribution among these stars would 
favor the thick disk having formed via the kincmatical heat- 
ing of an early thin disk or directly from merger debris. 

The Galactic metallicity distribution is best probed 
directly through spectroscopic surveys (Yoss, Neese & 
Hartkopf 1987; Allende-Prieto et al. 2006, hereafter AP06). 
However, an alternative method to spectroscopy is the use 
of field star UBV photometry (Majewski 1992; Gilmore & 
Wyse 1985; Karaali et al. 2003). UBV photometry allows 
the measurement of approximate abundances because the 
ultraviolet flux of a star is dramatically affected by the line- 
blanketing of the star's heavy metals - a phenomenon first 
described by Roman (1954) and Wildey et al. (1962). The 
measurement of ultraviolet excess - the amount of "extra" 
ultraviolet light produced by metal-poor stars in comparison 
to metal-rich stars of similar spectral type - is necessarily 
less precise than spectroscopy. However, it has the advantage 
of allowing the simultaneous measurement of many more 
stars out to the magnitude limit of a photometric survey. 

The photometric and spectroscopic studios referenced 
above have consistently shown vertical abundance gradients 
ranging from —0.75 to —0.08 dex kpc~^, with the gradi- 
ent slowly flattening at the high-2 regions dominated by 
the halo. However, these gradients do not necessarily reflect 
abundance gradients within any Galactic population. They 
could also be produced by the transitions from one mono- 
metallic population to the other. Majewski (1992), for ex- 
ample, argues from a photometric and astrometric survey of 
the North Galactic Pole that the thick disk has no vertical 
metallicity gradient and a non-Gaussian metallicity distri- 
bution. AP06, using a spectroscopic survey of kinematically- 
selected stars, demonstrate that neither the thick disk nor 
the halo have an intrinsic vertical abundance gradient. They 
also find that while the halo has a broad range of abun- 
dances, the thick disk has a compact metallicity distribution 
with a peak metalhcity of [Fe/H] ~ -0.7. Du et al. (2004), 
using an analysis of BATC photometry, posit that all three 
old stellar populations (old disk, thick disk and halo) lack a 
vertical abundance gradient. 

As this paper was in preparation, Ivezic et al. (2008, 
hereafter 108) published an extensive study of ugr photom- 
etry taken from the Sloan Digital Sky Survey. The tremen- 
dous number of stars available allows them to probe both 
broad and fine structure in the abundance distribution. 
While they confirm the general properties of the Galactic 
populations depicted in earlier surveys such as AP06, they 
do so with much more precision. Remarkably, their results 
would be consistent with the complete absence of a thick disk 
in the Milky Way or its reduction to the Metal- Weak Thick 
Disk described by Norris (1987). Karatas et al. (2008) have 
also recently analyzed ugriz photometry from the CFHTLS 
DR4. In contrast to other studies, they find a vortical abun- 
dance gradient in the halo. 

In this paper, we present an analysis of stars in the 
Kapteyn Selected Area SA 141, which is located near the 
South Galactic Pole and is ideal for probing vertical spa- 
tial, kinematical and abundance distributions. Trefzger, Pel 
& Gabi (1995) previously studied the photometric proper- 
ties of this field in an attempt to constraint the metallicity 



distribution. Using 112 stars over 1.9 square degrees to a 
depth of V — 14.8, Trefzger et al. found evidence for a dis- 
tinct thick disk population and gradients roughly consistent 
with previous spectroscopic studies. 

We now expand upon their result, presenting an anal- 
ysis of photoelectric UBVRI photometry in a similar field 
but to greater depth. The present study provides am inde- 
pendent check on the remarkable results of 108 as well as 
an exploration of issues affecting the derivation of metallic- 
itics from broadband UBV colours. Section 2 of this paper 
describes the observations, corrections for foreground extinc- 
tion and removal of QSOs contamination. Section 3 describes 
the method used to measure photometric abundances, ab- 
solute magnitudes and distances while Section 4 shows our 
analysis of the abundance distribution perpendicular to the 
Galactic disk and the comparison of our measures to those 
expected from the most recent models. 



2 OBSERVATION AND DATA REDUCTION 

2.1 Photometric Data 

The SA 141 field was observed as part of a larger starcounts 
survey, described in detail in S02. In brief, SA 141 was one 
of twelve fields observed in the UBVRI passbands with the 
wide-field Cassegrain CCD camera on the 1-m Swope tele- 
scope at Las Campanas Observatory between 1993 and 1998. 
Each field was observed in a uniform grid of pointings. In the 
case of SA 141, the U data cover approximately 1.3 square 
degrees to a mean depth of ?7 ~ 21. Transformation equa- 
tions to the standard system of Landolt (1992) were derived 
through matrix inversion techniques and each field was iter- 
atively transformed to this system. We applied the morpho- 
logical object classification techniques of S02 and identified 
1294 star-like objects with measurable U magnitudes. 

Since the publication of S02, we have slightly improved 
the data, generating astrometry from the IRAF TFINDER 
program and USNO SA2.0 catalog (Monet al. 1996). We 
also combined stars with multiple observations and removed 
slight photometric inconsistencies between the individual 
pointings using techniques described in Siegel et al. (2009). 

Photometry was de- reddened using the reddening maps 
and extinction coefficients of Schlegel, Finkbeiner & Davis 
(1998). Although several studies have indicated that the 
Schlegel maps over-estimated the reddening (see, e.g., Cam- 
bresy et al. 2005), the stars in the SA 141 field have minimal 
extinction (0.011 ^ E{B — V)^ 0.024), minimizing any im- 
pact on our study. 

2.2 QSO and galaxy contamination 

Deep, high-latitude fields can suffer significant contamina- 
tion from background galaxies. While S02's morphological 
techniques effectively remove most diffuse galaxies, compact 
galaxies can remain in the sample (Reid et al. 1996; S02). 
Fortunately, UBV photometry provides an excellent sec- 
ondary resource with which to remove compact background 
galaxies. 

Figs, la and c show the {U — B)o,{B — V)o and 
{V — R)o,{B — V)o colour-colour diagrams of SA 141. The 
sample shows a very clear and narrow stellar locus but has 
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Figure 1. Colour-colour diagrams for SA 141 field stars. Left panels (a) and (c) show (U — B)o,{B — V)o and (V ~ R)o,{B — V)o two 
colour diagrams of 1294 stars and star-like objects. Right panels (b) and (d) show {U — B)o, (B — V)o and {V — R}o,{B — V)o diagrams 
for 951 stars after possible galaxies and quasars are eliminated. 



significant contamination from compact galaxies and QSOs, 
manifested in the broader distribution overlying the narrow 
stellar locus. To remove these objects, we applied the selec- 
tion criterion of Chen et al. (2001), which removes objects 
with {u — g)o > 0.50. This transforms to a, U — B colour of 
—0.48, using the equations of Smith ct al. (2002). Wo further 
applied methods outlined in Fan ct al. (1999) and Karaali 
et al. (2003), which remove objects based on their location 
relative to the dominant stellar colour— colour locus. 

This cleaning reduced the sample from 1294 objects to 
951 stars (Figs, lb and d). The 343 objects removed from the 
sample is a far greater than the level of QSO contamination 
expected (sec, e.g., Richards ct al. 2001) and it is likely that 
our photometric locus method is removing outlier stars that 
have compromised measures or unusual properties. However, 
the objects removed from our analysis lie well outside the 



boundaries of the iso-metallicity lines discussed in Sect. 3.3 
and their absence does not affect our analysis. 

3 ANALYSIS TECHNIQUES 
3.1 Stellar population types 

The (Vo, {B — V)o) colour-magnitude diagram for the 951 
stars in the SA141 sample is plotted in Figure 2. We have 
indicated the cut-off colour indices for thin and thick disks 
and halo populations as described by Chen et al. (2001, fig- 
ure 6). These cutoffs occur at {g — r)o colours of 0.33 and 
0.20, respectively, which wc convert to {B — V)o colours 
of 0.54 and 0.40 using the transformation of Smith et al. 
(2002). Our photometric distribution is roughly consistent 
with Chen, with the thick disk dominating the starcounts 



4 Siegel et al. 



12 



13 



14 



15 - 



r.i 



• • • 



Thick-disk cut 



16 - 



> 17 



, ' . • . . .• 

• • J— i • • • ' 
.'• • • 



18 - Halo cut-off 



20 - 



21 - 



22 




-0.5 



0.5 



1.5 



(B-V)3 



Figure 2. The {B—V)o — Vo colour-magnitude diagram of 951 
stars in SA 141 field. Dashed lines show the cut-offs of (B — V)o 
for the thick disk and halo derived by Chen et al. (2001). 



to about Vb ~ 19 and the halo dominating at fainter mag- 
nitudes. 



3.2 Model Comparison 

Untanghng stellar populations from raw colour-magnitude 
diagrams is a difficult process at best and frequently prone 
to degenerate solutions, especially when based upon obser- 
vations in a single direction. Sophisticated models, which 
incorporate up-to-date spatial and abundance distributions, 
luminosity functions and isochrones, can be used to make 
Monte Carlo comparisons with which to constrain the prop- 
erties of the underlying stellar populations. The most readily 
available and recently updated model is the Besancon model 
of stellar population synthesis (Robin et al. 2003, hereafter 
which can produce synthetic photometric catalogs 
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based on an updated Galactic structure model and an input 
set of magnitude limits and error functions. 

Figure 3 compares the Besancon simulation to the SA 
141 photometry. For the simulation, we used the updated 
default structural parameters given by the Besancon simu- 
lator. We simulated errors using exponential distributions fit 
to the mean error locus of the data (Figure 4). No a priori 
color limits were included and the magnitude limit was set 
to the mean imaging limit of the SA141 data {V = 21.2; 
B = 20.9; U = 21.0). The imaging limit was defined in 
S02 as the faintest magnitude at which stars and galaxies 
could be distinguished based on DAOPHOT morphological 
parameters. 

Both the real and simulated data fall along an arced lo- 
cus in colour-colour space, with the location of each star be- 
ing a function of its effective temperature and stellar metal- 
licity (as explored in detail in Sect. 3.3). The numerous blue 
stars that create a clump in upper left corner of the colour- 
colour diagram would correspond to the blue edge depicted 
in Fig. 2 and shown in the lower panels of Fig. 3. The blue 
edge is comprised of intrinsically bright early-type stars. 
These stars are faint in the CMD and therefore several kpc 
distant in regions dominated by the metal-poor halo. The 
long trail of stars to redward is the increasing contribution 
of the intrinsically fainter (and therefore nearer) late-type 
stars, which would be in regions of the Galaxy dominated 
by the thick and thin disks. 

While the basic features are similar, a number of dis- 
crepancies are seen. The simulated blue edge is much more 
populous than the observed blue edge. The simulated blue 
edge also has a much tighter photometric distribution than 
the real blue edge in both colour-magnitude and colour- 
colour space. Since the Besancon model accounts for photo- 
metric error and reproduces the redder features of the CMD 
quite well, it is likely that this discrepancy is intrinsic to the 
model. 

The Besancon model uses a halo comprised of a sin- 
gle stellar population which has a simple power law spatial 
distribution, a mean abundance of [Fe/H] ~ —1.78 and 
a metallicity dispersion of 0.5 dex. However, the numerous 
studies referenced in Sect. 1 indicate that simple models 
have difficulty reproducing the halo spatial distribution and 
demonstrate that the halo is partly, perhaps mostly, com- 
prised of streams of stars stripped from objects during the 
hierarchical formation of the Milky Way. The stark differ- 
ence in uniformity between the real and simulated blue edge 
argues against simple single-population models of the deep 
halo. 

We also note that the simulated galaxy has a slightly 
different overall colour-colour distribution, with a gentler 
slope and a bluer cutoff in the {U — B)o colour distribu- 
tion. This may not indicate a model deficiency as much as 
a mismatch between the real and simulated filters. Our ob- 
servational data are calibrated to the standards of Landolt 
(1992) while the Besancon simulation is tied to the model 
atmospheres of Lejeune, Cuisinier & Buser (1997, 1998), 
which itself is tied to the observational plane by the calibra- 
tion of Schmidt-Kaler (1982). It is unsurprising that some 
discrepancy would arise between two different calibrations. 
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Figure 3. A comparison of the SA141 piiotometry (left panels) to the Besancon simulation (right panels). Note the much tighter blue 
edge population in the simulation. 



3.3 Metal abundances and absolute magnitudes 

The metal content of main sequence stars can be mea- 
sured by the amount of line-blanketing absorption affect- 
ing their ultraviolet photometry. Low metal abundance pro- 
duces weak atmospheric absorption lines and increased ul- 
traviolet emission, or ultraviolet excess (UVX; Wildey et al. 
1962). 

Sandagc (1969) detailed a technique for using the UVX 
- as measured through UBV photometry - to measure ap- 
proximate abundances. The method determines the two- 
colour location of a star with reference to iso-metallicity 
ridgelines. The UVX measure, 5{U — B), is the difference 
between each star's {U — B)o colour and that of a metal- 
rich star of identical {B — V)o, with the reference popu- 
lation defined in this case by the Hyades sequence, which 
Cayrel, Cayrel de Strobel & Campbell (1985), Boesgaard & 
Priel (1990) and Taylor (1998) estimate has an abundance 
between [Fe/H]^+O.W and [Fe/ H]=+0.13. We adopt an 
abundance of -1-0.13 dex for our analysis. 

Figure 5 shows the de-reddened two-colour, {U-B)o- 
{B-V)o of the SA 141 stars that fall within the colour range 
of the Karatas & Schuster (2006) UVX calibration. The solid 
lines indicate the metal content, ranging for the reddest line, 
which is a metal-rich line from the Hyades sequence, to the 
bluest line, which indicates the maximum underabundance 
envelope defined by Sandage (1969). 

Inspection of Fig. 5 shows that the stars at bluer 
{B — V)o colours tend to lie along the more metal-poor 
lines while redder stars tend to lie along the more metal- 
rich lines. This is in keeping with expectations of standard 
Galactic structure models. The redder stars are intrinsically 
faint and nearby, in regions dominated by the metal-rich 



thin disk. The bluer stars arc intrinsically bright and distant, 
in regions dominated by the metal-poor halo. A significant 
number of faint blue halo stars lie beyond the most metal- 
poor lino, which corresponds to [Fe/ H] =—2.23 in the cali- 
bration of Karatas & Schuster (2006). This could represent 
the low-metallicity tail of the halo but would also represent 
photometric scatter and evolved stars that have colours out- 
side the main sequence locus. A smaller number of nearby 
red disk stars lie rodward of the Hyades line, possibly rep- 
resenting the high-metallicity tail of the disk. Excluding the 
outlying stars from analysis can lead to a systematic bias, 
which we address in Sect. 4.2. 

To derive [Fe/H] abundances, we first measure 6{U — 
B), the colour difference between each star and the maxi- 
mum metallicity line. This measure is then transformed to 
6o.6, the UVX measure normalized to a {B—V)o colour of 0.6 
(where S{U — B) reaches its peak value) along the indicated 
iso-metallicity ridgelines. Eight constant-abundance lines of 
Sandage (1969), drawn in Fig. 5, have been used to trans- 
form the 5{U-B) measures of our stellar sample into nor- 
malized 6o.6 measures. Constant-metallicity lines are scaled 
as a fraction of the maximum line-blanketing "M" value in 
increments of 0.125. The lines are taken from Table lA of 
Sandage with adjustments as detailed in equations (2)— (4) 
of Karatas & Schuster (2006). 

We convert the So.e measures to approximate abun- 
dances using the Karatas & Schuster (2006) relation: 

[Fe/H] = -h0.13(±0.04) - 4.84(±0.60)5o.6 - 

7.93(±2.24)5^.6 (1) 

This relation is valid for the ranges of -1-0.41 < [Fe/H] ^ 
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Figure 4. Error measurements for stars from the SA141 field. The error trends show multiple loci owing to the varying quality of the 
data. The solid line represents the error function of the simulated data. 



-2.23, -0.04 < 5o.e> < +0.29 and 0.3 < (B - V)o 0.9. 368 
of the SA 141 stars fall within this parameters space. 

Absolute magnitudes of program stars are then esti- 
mated with the following relation of Karatas & Schuster 
(2006): 

Mv = 2.77 + 4.38(B-V)o - 48.985o.6 
+ 219.75(B-V)o Sle - 198.93(B-V)o 5o.6 

+ 11.04(B-V)2 So.6 (2) 

Abundance uncertainties are calculated by propagating 
the photometric uncertainties through the relevant equa- 
tions. The abundance uncertainties are large, starting at a 
mean 0.2—0.3 dex for nearby {z < 4 kpc) stars, reaching 
0.5 dex at 8 kpc and 0.7—0.8 dex at 15 kpc. This indicates 
that the abundance distribution of the thin and thick disks 
should be well-constrained. However, the deep halo, where 
the number of stars is low and the photometric uncertain- 
ties high, may be poorly constrained, an issue we address in 
Sect. 4.2. 

Our measurement of stellar abundances, absolute mag- 
nitudes and distances is predicated on the assumption that 
these stars are all main sequence staxs. This is not the case 
and some degree of contamination from evolved stars is ex- 
pected. The evolved stars will be more distant than the main 
sequence stars and consequently more metal-poor, given the 
abundance distribution identified by previous investigators 
and constrained in Sect. 4. This may cause analysis of the 
nearby sample to skew metal-poor. 

Fortunately, the Bosancon simulated galaxy retains in- 
formation on luminosity clciss. To evaluate the effect of this 
on our derived abundance distribution, we ran the simu- 
lated galaxy through our analysis pipeline and examined 



UVX-based metallicity measures as a function of midplane 
height, both with evolved stars removed from the sample 
and with evolved stars retained and evaluated as though 
they were main sequence stars. We found that retaining the 
evolved stars caused a small metal-poor skew in the derived 
metallicity that had a maximum of —0.1 dex but was goner- 
ally around a few 0.01 dex. The reason the effect is so small 
is that while the evolved stars comprise 10-20% of the star- 
counts at any magnitude, a significant number lie outside the 
maximum and minimum UVX boundaries shown in Fig. 5 
and arc rejected on that basis. Indeed, it is likely that some 
of the outliers in Fig. 5 are, in fact, evolved stars. Addition- 
ally, the relatively blue cutoff of our abundance and absolute 
magnitude calibrations removes the redder stars where the 
luminosity difference between giant and dwarf is larger and 
the contamination is therefore from far more distant and 
more metal-poor stars. We have chosen not to remove the 
evolved star bias as it is smaller than the abundance uncer- 
tainties. 

4 METALLICITY DISTRIBUTION 
4.1 Raw Abundance Gradients 

The trend of mean metal abundance < Fe/H] > against 
< 2 > height above the Galactic plane is shown in Figure 6 
and detailed in Table 1. Mean metal abundances are derived 
from maximum likelihood Gaussian fits for star numbers 
and z intervals given in Table 1. As can be seen in Figure 
6, at distance of z < 4 kpc - where the S02 and Robin et 
al. models indicate that the thin and thick disks dominate 
- there is a vertical abundance gradient of approximately 
d[Fe/H]/dz = —0.15 dex kpc~^. At distance beyond 5 kpc 
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Figure 5. (U-B)o-{B—V)o two-colour diagram for stars in SA 
141 field. Standard relations for the Hyades and for the metal- 
poor lines of (!7-S)o.25M, (C^--B)o.50M , {U-B)o,y5M, and [U- 
B)i^ (metal-free) from Sandage (1969) are drawn as solid lines. 
Intermediate metal-poor lines of ((7-i3)o.i26M) {U-B)q 3ysm, 
({7-B)o.625M , {C^-B)o.875M are drawn as dashed lines. The deci- 
mal fractions indicate the proportion of the maximum ultraviolet 
excess, which corresponds to zero metallicity. Open points are 
stars that fall outside of our constrained abundance range. 

- where the halo begins to dominate - the metallicity mear 
sures are less precise but are consistent with a flattening 
of the abundance gradient (0.03—0.06 dex kpc^^). At large 
distances (> 12 kpc), the abundance seems to rise slightly. 
However, the data at these faint magnitudes suffers from low 
completeness levels and it is likely that this reflects a bias 
near the magnitude limit or contamination from evolved hor- 
izontal branch (HB) stars (Sections 4.2 and 4.3). 

4.2 Debiasing and Deconvolving the Abundance 
Distribution 

The raw mean abundances shown in Figure 6 can not be 
taken at face value due to two complicating factors. First, 
the vertical abundance distribution is a convolution of the 
intrinsic abundance and density distributions of the compo- 
nent populations. Second, the method itself induces bias. 
Our metallicity calibration only extends from [Fe/H] = 
-1-0.41 to [Fe/H] = -2.23. This cuts off the low-metallicity 
tail of the halo and the high-metallicity tail of the disk - a 
process shown visually in Figure 5. Removing these outlying 
stars results in a systematic bias. 

The dramatic effect of this bias can be demonstrated by 
applying our UVX analysis pipeline to the Besancon sim- 
ulated data. Figure 7 compares the simulated abundances 



Table 1. Mean metal abundances for z intervals (Column 1). 
Column 2 shows the mean < 2 > of the stars in each bin while 
Columns 3 and 4 show the mean < [Fe/H[> and uncertainty of 
the mean, respectively. Column 5 indicates the calculated abun- 
dance dispersion of the stars while Column 6 indicates the number 
of stars in each bin. 



z range 


<z> < [Fe/H] > 


<^<[Fe/H]> 




N 


[0,1] 


0.64 -0.50 


±0.07 


0.56 


75 


(1,2] 


1.48 -0.61 


±0.06 


0.55 


116 


(2,3] 


2.47 -0.79 


±0.12 


0.59 


43 


(3,4] 


3.53 -0.93 


±0.11 


0.28 


31 


(4,5] 


4.52 -0.65 


±0.17 


0.48 


23 


(5,6] 


5.38 -0.96 


±0.21 


0.60 


18 


(6,7] 


6.52 -0.95 


±0.23 


0.31 


12 


(7, 8] 


7.34 -0.55 


±0.24 




9 


(8,9] 


8.52 -1.23 


±0.35 




8 


(9, 10] 


9.38 -1.13 


±0.29 




9 


(10, 11] 


10.54 -0.99 


±0.33 




7 


(11, 13] 


12.31 -0.27 


±0.25 




9 


(13, 15] 


14.02 -0.43 


±0.20 




11 


[15,21] 


16.73 -0.50 


±0.39 




7 
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z (kpc) 

Figure 6. Relation between < [Fe/H\ > and < z > for the 
range of ^ z ^ 20. There is a slight metallicity gradient 
d[Fe/H]/dz = —0.15 dex kpc~^ for 2 < 4 kpc, where the thin and 
thick disks dominate. Error bars in the top panel reflect the for- 
mal uncertainty in the mean calculated by maximum likelihood 
from the input photometric errors. Error bars in the bottom panel 
reflect the measured intrinsic abundance dispersion. 

with the UVX-derived abundances. While the two scales fol- 
low each other closely at the metal-rich end, the trend veers 
at low abundances. The reason is that the metal-poor halo 
stars have a large abundance dispersion - both intrinsic (the 
Besancon halo has a 0.5 dex dispersion) and observational 
(0.5—0.8 dex from the UBV photometric uncertainties prop- 
agated through the UVX equations). The low-abundance 
cutoff removes the (intrinsic plus photometric scatter) low- 
metallicity tail of this distribution, while leaving the high- 
metallicity tail intact. In the faintest reaches of our halo 
sample, this shifts the derived mean abundance by as much 
as 0.4 dex. 

Figure 8 contrasts the underlying intrinsic mean abun- 
dance distribution against the UVX-derived mean abun- 
dance distribution. The difference between the two is stark 
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Figure 7. A demonstration of the bias induced by our UVX 

analysis. The points show the input abundances from the Besan- 
con simulation contrasted against those derived from the UVX 
pipeline and simulated UBV photometry. The line represents 
unity. The sudden shift of the low abundances away from unity 
is not the result of miscalculation or error but the increase in dis- 
persion combined with the removal of the metal-poor measures 
whose UVX properties lie beyond the last ridgeline in Figure 5. 




Figure 8. A demonstration of the metallicity bias induced by 
UVX analysis. The circles represent the intrinsic mean abun- 
dance distribution of the Besancon simulated data. The square 
points represent the derived mean abundance measures from the 
simulated UBV photometry. In both cases, the uncertainties are 
calculated from maximum likelihood (with an assumed uncer- 
tainty of a-p^/fj = 0.2 for the intrinsic abundances). Because 
the UVX method has large uncertainties at faint magnitudes and 
only extends to [Fe/H] = —2.23, a significant metal-rich bias is 
induced when the stars scattered to low abundances are removed, 
as demonstrated in Figure 7. 



and is the result of the abundance cut-off. A similar bias 
occurs for the metal-rich thin disk, but is not as dramatic 
owing to the smaller uncertainties for the bright disk stars 
and the smaller (0.2 dex) intrinsic dispersion of thin disk 
abundance. Note also that the error bars in Figure 8 are 
smaller than those of Figure 6, owing to the larger number 
of faint halo stars in the real data. 

Untangling the halo metallicity distribution under these 
circumstances is beyond the scope of this paper -and may in- 



deed be impossible if the halo metallicity distribution is not 
a simple single population but, as seems likely, comprised of 
multiple components and/or stream of stars. However, since 
(a) the SA 141 data and UVX calibration are on the same 
photometric system (both are tied to the Landolt standards) 
and; (b) the nature of the low-metallicity bias is known, it 
is possible to construct toy models for comparison. 

We constructed synthetic metallicity distributions by 
assuming the structural parameters detailed in S02, metal- 
licity distributions from AP06 and R03 and a distance- 
uncertainty relation taken directly from the SA 141 data. 
Both formulations posit a halo and thick disk with no metal- 
licity gradient and abundance dispersions of 0.3 and 0.5 dex, 
respectively. The AP06 mean metallicities arc —0.68 and 
— 1.4 dex for thick disk and halo, respectively. The corre- 
sponding R03 abundances are —0.78 and —1.78 dex. For 
the thin disk, we took the Besancon metallicity distribu- 
tion, which is a combination of seven thin disks in which 
metallicity is dependent on age. The mean metallicity of the 
disks varies with age. However, the overall weighted mean 
abundance of the seven disks is [Fe/H] = —0.08 and the 
weighted weighted Icr is 0.2 dex. We applied the low- and 
high- abundance cutoffs that are applied to the real data to 
replicate the bias. 

Figure 9 shows the result of these comparisons. What 
is striking about the distributions is that the photometric 
properties of the Galactic populations show a rough distri- 
bution consistent with the models, but overlayed with sev- 
eral metal-rich outlying data points (at 4.5, 7.5 and beyond 
12 kpc). S02 (and many studies cited therein and in Sect. 1) 
have suggested that the halo field star distribution has two 
components. The first is a uniform fiattcncd metal-poor dis- 
tribution; the second an irregular mass of star streams cre- 
ated during the hierarchical formation of the Milky Way 
(although Carollo et al. 2007, using the SDSS, argue for 
two distinct halo structures, only distinguishable beyond 10 
kpc) . If the halo did have a uniform structure overlayed with 
streams, the abundance distribution in a pencil-beam survey 
might look like Figure 9, with a uniform distribution occa- 
sionally interrupted by more metal-rich populations. Inspec- 
tion of Figure 6 shows a hint of this at the outlying points, 
where the individual points seems to show a vague bimodal- 
ity, although the scatter it too large to be certain. 

Deep proper- motion and/or radial velocity surveys of 
SA 141 would be needed to test this suggestion. It remains 
possible that those data points merely reflect some inad- 
equacy of the method. Nevertheless, if verified, it would 
demonstrate that the stellar streams identified so clearly 
in 108 can also be identified in smaller less-precise archival 
UBV data sets by their departure from the underlying uni- 
form distribution. 

If we assume that the metal-rich points are outliers 
above the general distribution, then the literature models 
successfully reproduce the observed abundance distribution. 
The Besancon model appears to provide a slightly better de- 
scription of the thick disk while the AP06 model provides a 
slightly better description of the halo. The halo is not well- 
constrained enough to address a possible abundance gradi- 
ent. A slight thick disk abundance gradient would bring the 
AP06 thick disk in line with the UVX results, albeit only 
because it would lower the mean abundance in the narrow 
range of z heights dominated by the thick disk. 
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Figure 9. Comparison of tlie UVX-based abundances to tliose 
predicted by tlie Besancon and AP06 distributions (solid and dot- 
ted lines) as well as both distributions with a —0.05 dex/kpc 
abundance gradient in the thick disk (dashed and dot-dashed 
lines). In this figure, the model predictions do not correspond 
to the actual metallicity distribution, but to the measured abun- 
dance distribution, taking into account the bias induced by the 
UVX methodology. 

4.3 Young Stars in the Halo? 

Figure 10 shows the [Fe/H], {B-V)o diagram overlayed with 
Yonsei-Yale (Y^) isochrones of Yi, Kim &Demarque (2003). 
The turn-off loci from the isochrones are overlayed for ages of 
6, 8, 10, 12, and 14 Gyr, assuming [a/Fe] = +0.30. While 
the vast majority of the stars in SA141 have photometric 
properties consistent with ages greater than 12 Gyr, a hand- 
ful of metal-rich stars are bluer than the old turnoff. These 
stars are all beyond 6 kpc and would have implied ages less 
than 10 Gyr. A young metal-rich population of stars is un- 
expected given previous spectroscopic and photometric sur- 
veys of the halo (see CaroUo et al. 2007). 

None of the stars comprising this young population have 
particularly large photometric errors. Nor are they concen- 
trated in one part of the SA 141 field, eliminating the pos- 
sibility that they are compromised by bad pixels or nearby 
saturated stars. Although the abundance uncertainties are 
large for the fainter stars (0.5-0.8 dex) and they suffer from 
the aforementioned metal-poor bias, the mean abundance is 
clearly metal-rich (-0.2 with a dispersion of 0.7). 

These stars are unlikely to represent the halo streams 
to which we ascribe the outlier points in Figures 6 and 9. 
Halo stream stars are generally thought to be metal-poor - 
as reflected in the blue, metal-poor (BMP) stars described 
in Preston, Beers & Shectman (1994), Carney et al. (1996) 
and Unavane, Wyse & Gilmore (1996). While some streams, 
most notably those of the merging Sagittarius dSph, are 
known to be slightly more metal-rich than the halo (Sbor- 
done et al. 2004; Martinez-Delgado et al. 2005; Chou et al. 
2007), our distant SA 141 stars are not clumped in distance, 
as would be expected for a tidal stream. 

It's possible that these stars are blue straggler stars 
(Carney, Latham & Laird 2005) although there is no rea- 
son to believe that blue stragglers would be so much richer 
than the canonical halo. However, the most likely explana- 
tion is that these are contaminating HB stars. Figure 11 
shows the SA 141 UBV diagram with the young stars of 
Figure 10 marked as large circles. To contrast their loca- 
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Figure 10. (a) The [Fe/H], {B-V)o plot for the 368 stars in 
SA 141 field. The turn-off loci from the isochrones of Yi et al. 
(2003) are plotted for the ages of 6, 8, 10, 12, and 14 Gyr, with 
[a/Fe] = -1-0.30. The lines indicate the young metal-poor stars 
surveyed by Unavane et al. 

tion against that of HB stars, we have generated a synthetic 
HB from the Dartmouth Stellar Evolution Database (Dot- 
ter et al. 2007), assuming an abundance of [Fe/H] — —1.5, 
[a/Fe]— +0.2, an age of 12 Gyr and mass loss of 0.05 ± 
0.05 McQ. It can be seen that the HB somewhat overlaps 
the metal-rich end of MS in colour-colour space. If the dis- 
tant blue metal-poor "young" stars were on the HB, rather 
than the MS, they would be spread over distances of 25—150 
kpc, a distance range dominated by the extremely metal- 
poor {[Fe/H] ~ —2.2) outer halo described by CaroUo et 
al. (2007). We find this the most likely explanation for these 
outlying stars. Spectroscopy or proper motions of these stars 
would be able to determine their nature. 



5 CONCLUSIONS 

We have performed an analysis of the ultraviolet excesses 
of stars in the SA141 field. We find that photometric un- 
certainties propagate very quickly into large (0.5—0.8 dex) 
abundance uncertainties. These uncertainties can cause sys- 
tematic biases at the low-metallicity end of the UVX cali- 
bration. These biases must be properly accounted for when 
probing the ensemble properties of the field stars. Specifi- 
cally, the limitation of the method to stars more metal rich 
than [Fe/H] = —2.23), when combined with the intrinsic 
abundance dispersion of the halo and the necessarily large 
UVX-abundance uncertainties at faint magnitudes, induces 
a systematic bias in our analysis of the population abun- 
dance. This bias can be modeled and corrected if assump- 
tions are made about the underlying metallicity distribution. 

^ The location of the SHB is only sensitive to the input parame- 
ters at colours bluer than the colour range over which we measure 
UVX abundances. 
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Figure 11. A analysis of the unusual halo stars. The dots show 
the U BV photometry of the SA141 field stars. The large circles 
are the outlying apparently young stars from Fig. 9. The triangles 
are a synthetic horizontal branch. Note the overlap of the HB and 
MS sequences where the outlying stars are found. 



The inferred abundance distribution is roughly consis- 
tent with the metallicity distribution depicted in 108 as well 
as the spectroscopic study of AP06 and the simulations of 
R03. However, there are a number of outlying data points 
that could indicate contamination of the SA 141 field by 
more metal-rich halo streams. 

We identify a trace population of apparently young 
metal-rich halo stars. However, these are not associated with 
any of the outlier data points and we find it is most likely 
that these objects are halo HB stars that have photometric 
properties similar to metal-rich MS stars. 

To exploit future UBV data-bases, it will be necessary 
to refine the calibration of the UBV ridgelines based on 
up-date photoelectric data calibrated to the Landolt system 
as well as to probe the UBV properties of very metal-poor 
stars to extend the metallicity range of the calibration. This 
would allow the abundance distribution of the distant halo 
to be explored with more precision and confidence to test 
models of Galactic formation and unravel the fine structure 
of the Galactic halo and thick disk. 
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